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Abstract 

We show that charged black holes in Anti-dc Sitter spacctimc can undergo a third 
order phase transition at a critical temperature in the presence of charged fcrmions. In 
the low temperature phase, a fraction of the charge is carried by a fermion fluid located 
a finite distance from the black hole. In the zero temperature limit the black hole is 
no longer present and all charge is sourced by the fermions. The solutions exhibit the 
low temperature entropy density scaling s ~ T 2 / z anticipated from the emergent IR 
criticality of recently discussed electron stars. 



Recent work has argued that 'electron stars', a planar fluid of charged fermions in Anti- 
de Sitter space in equilibrium under gravitational and electromagnetic forces, provide a 
compelling holographic framework in which to study metallic quantum criticality [USE]- 
Two key features of electron star solutions are emergent criticality at low energy with a 
finite dynamical critical exponent [JJ [2] and the presence of a 'smeared' Fermi surface [3]. 
Metallic criticality is difficult to study with conventional field theoretic techniques as the 
many gapless excitations of the Fermi surface cannot be integrated out and must be included 
in a strongly interacting IR fixed point [I]. 

The zero temperature electron stars studied in [2] are also interesting as charged 'soli- 
tonic' gravitational configurations without a black hole horizon. Their existence is aided by 
the gravitational well of the asymptotically Anti-de Sitter spacetime (AdS). In this sense 
they are the fermionic analogues of zero temperature holographic superconductors El [7] . 
At sufficiently high temperatures in the holographically dual nonzero density field theory 
we might expect the electron star to collapse and form a Reissner-Nordstrom black hole. 
This would be analogous to the fact that zero temperature stars in AdS with spherical, 
rather than planar, symmetry undergo a first order transition to a black hole as a function 
of the energy of the star in units of the size of the spatial sphere at the boundary [U [9]. 
The primary result of this letter is that electron stars undergo a third order transition to a 
charged black hole above a critical temperature determined by the chemical potential and 
mass of the fermions. Lowering the temperature through the transition, the birth of the 
electron star, thereby corresponds to a continuous phase transition in the nonzero density 
dual field theory characterised by the appearance of a smeared Fermi surface. 

This letter will study the dynamics of a 3+1 dimensional ideal fluid of charged relativistic 
fermions coupled to electromagnetism and gravity with a negative cosmological constant. 
The Lagrangian may be written 



This is the Schutz form of the action for a gravitating ideal fluid in terms of the fluid 
pressure p |10| , generalised to allow the fluid to be charged [2] . The local chemical potential 
[i = \dcj) + ad(3 + Ods + A\, where {(ft, a, /?} are fluid potential variables, s the local entropy 
density and 6 the thermasy. This action leads to the expected ideal fluid equations of motion 
as described in e.g. [T01 12"]. 

The fluid Lagrangian Q is a coarse grained description of the fermions in which the 
fermion physics is subsumed into locally defined thermodynamic quantities. In the context 
of astrophysics this could be called the Tolman-Oppenheimer-Volkoff description [Til 02] 
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while in condensed matter physics it is known as the Thomas- Fermi approximation [1 3|. HH] . 
In our context, the fluid Lagrangian ([I]) is a correct description of the system when [2] 

e 2 ~ y < 1 . (2) 

Without loss of generality we will take the fermions to have unit charge and mass m. 
For detailed discussions of the connection between microscopic and fluid descriptions of 
gravitating fermions, see e.g. [T5j 19], 

The spacetime metric and Maxwell field take the general form 

ds 2 = L 2 {-fdt 2 + gdr 2 + 1 (dx 2 + dy 2 )^j , A = ^hdt . (3) 

For the fermion fluid, a crucial role is played by the local chemical potential 

A t e h 

^ =A * = WJ = -.7J' <4) 

This chemical potential determines the local thermodynamic quantities of the fermion fluid. 
Before specifying these, however, we need to discuss how a nonzero temperature is going 
to appear in our equations. Placing the dual field theory in equilibrium at a nonzero 
temperature can be described by compactifying the Euclidean time direction. This implies 
that the Euclidean time direction in the bulk must also be compactified. Two consequences 
of a periodic bulk Euclidean time circle are that for regularity of the spacetime, recall that 
we are in planar coordinates, we must have a finite size black hole horizon in the interior and 
secondly that the local fermion fluid equation of state is at finite temperature. Physically 
this describes a black hole surrounded by a fermion fluid in thermal equilibrium with the 
Hawking radiation. In our bulk classical limit the effects of Hawking radiation should be 
negligible, while the black hole remains present. Therefore, we should expect that in this 
limit we can treat the fermions as a zero temperature fluid in a black hole background. We 
will make this statement precisely towards the end of the paper. 

The upshot of the previous paragraph is that we can use the same zero temperature 
equations for the Einstein-Maxwell-fluid system as for zero temperature electron stars, but 
with boundary conditions ensuring that the spacetime has a finite temperature horizon. 
This means that the energy density, charge density and pressure of the fermion fluid are 
determined by the local potential Q via their zero temperature equation of state 

h h 

p = /3 e 2 \/e 2 — rh 2 de , & = fi I e\J e 2 — rh 2 de , ~P = P t=<t . (5) 

J fh J m V J 

Here we have introduced dimensionless hatted quantities 



Throughout we will use hats to denote bulk or field theory quantities that have been rescaled 
by factors of {e, L, k} in such a way that no such factors appear in the equations. The 
constant (3 is a dimensionless order one number that counts the microscopic degrees of 
freedom of the fermion. The important fact about $ is that it measures the ratio of the 
Maxwell and Newton couplings and should also be order one if we are in the regime ([2]). 
Furthermore, the equations of motion following from the action ([I]) imply the same equations 
for the components of the fields as in [2] 

1 ( f 9' 4\ aha 

+ - = 0, (7) 



r \f 9 rj 
f h' 2 1 

jj- y+9(Z + P)-- 2 = 0, (8) 



V7 V 2 

From ([5]) it is immediate that a nonvanishing density of fermions at a particular radius 
requires that the mass be lower than the local chemical potential m < ^i oc ., or m < hj\f~j '. 
In the Reissner-Nordstrom black hole solution 

1/1 fi 2 \ r (j? r 2 1 i ~ ( r 



f = ~2~ Hr + ^- — + ^-2-' 9 = — f , h = p, 1-— , (10) 

one can see that h/y/f is bounded with a maximum in between the horizon at r = r + 
and boundary at r = 0. Thus if m is too large, the fluid density will vanish everywhere, 



p = a = p = 0, and Reissner-Nordstrom (10) will be the exact solution. As the mass is 



lowered it will eventually equal the local chemical potential at a critical radius 

= 0. (11) 



h 
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d h 

m , = 

dry/f 



It is simple to solve these two equations. The two equations determine the critical radius 
at which the star is born, dimensionlessly expressed as r c /r+, and the critical temperature 
of the black hole over the rescaled chemical potential Tc/fi- Recall that the temperature of 
a Reissner-Nordstrom black hole is 



1 

Aire 



( 1 
dr 



(12) 



At this point c = 1, we have included this factor for later convenience. 

In figure[T]we plot the critical radius and critical temperature as a function of the fermion 
mass. We see that the electron star forms at increasing low temperatures, and closer to the 
black hole horizon, as the fermion mass is increased. When the mass is increased beyond 
m = 1 the star is never formed. This is consistent with the observation in [2] that rh < 1 is 
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Figure 1: Critical temperature (left plot) and radius (right plot) at which the electron star 
is born, as a function of the fermion mass. These quantities do not depend on f3. 



necessary for the zero temperature electron star to exist. The critical temperatures and radii 
plotted in figure [I] follow from the Reissner-Nordstrom black hole geometry and therefore 
do not depend on the parameter /3 appearing in the fluid equation of state (ph. In fact, the 



conditions (11) simply characterise the radius at which a charged point particle can remain 



stationary, with gravitational and electromagnetic forces balancing. 

Cooling below the critical temperature of figure[I]we will find two solutions {r\, rg} to the 
equation /Lii oc . = m defining the boundary of the star. Thus the electron star broadens into 
a thin shell. Once broadened, the equation of state of the fluid as well as the backreaction 
of the fluid onto the geometry will be important. The solution at T < Tq has three 
components: 



1. Inner region, r > T2- The solution has the Reissner-Nordstrom form (10), but with 
fi — > pLQ, not directly related to the chemical potential of the dual field theory. This 



region determines the temperature through (12). 



2. Intermediate region, r% > t > T\. Here we must solve the equations ([7]) to ^ 
with p, a,p nonzero. The quantities {/, g, h, h'} must be matched onto the inner and 
exterior regions at and r%. 

3. Exterior region, r\ > r. The solution again has the Reissner-Nordstrom form, al- 
though now it should be written as 

c 2 

r 4f 



f 



h = c ( /* — rQ 



(13) 



The factor c, determining the normalisation of time, must be included in the definition 



of the temperature (12). Note the definitions are slightly modified with respect to [2]. 
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In performing numerics one may without loss of generality set r+ = 1. The solutions 
are then parametrised by £lq in the inner region. Via matching at r<i and r\, this initial 
condition then integrates forward to determine the values of the physical dual field theory 
quantities {c, M , Q, /t}. The solutions can then be labelled by the physical dimensionless 
ratio T/fi or alternatively T/Tq. Performing the numerics is simple and Figure [2] shows 
two examples of how the density of fermions a builds up and expands to widening shells as 
the temperature is lowered below the critical value. In reading these plots one should keep 




Figure 2: Radial density profiles of electron stars as a function of temperature. The curves 
show five temperatures between 0.07 Tc and Tq, with fi held fixed. Both plots have (3 = 10. 
The left plots has mass rh = 0.7 while the right plot has rh = 0.1. 



in mind that the total charge in the fermion fluid will be [2] 

Qformion = / —&(s)ds. (14) 

In particular, it is necessary to include the spatial volume element y/g/r 2 in summing up 
the contributions from different radii to the total charge in the fermion fluid. 

The total charge density of the boundary field theory, read off from the exterior space- 



time solution (13), may be written 

Q = QBH + Qfermion , (15) 

where the charge density carried by the black hole is read off from the interior solution as 

Qbh = ^ • (16) 
r+ 

From the bulk perspective, the presence of the fermionic charge density is the defining 
characteristic of the electron star. This suggests examining the temperature dependence of 
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the fraction of charge carried by the fermionic fluid, (Q — Qbr)/Q, as a sort of bulk order 
parameter. The ratio will be zero for temperatures above the critical temperature and, 
based on the observation in [2] that at zero temperature all the charge is carried by fermions, 
should tend to unity at low temperatures. Figure [3] shows precisely this phenomenon. For 
both the low mass and heavier fermions, the charge is transferred from the interior black 
hole to the fermion fluid more slowly as the temperature is lowered than for intermediate 
mass fermions. 
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Q 
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Figure 3: The fraction of charge carried by the fermion fluid as a function of tem- 
perature. All curves have /3 = 10 while from blue to red (solid to dotted) m = 
{0.01,0.07,0.3,0.55,0.75}. In quoting T/T c , the chemical potential is kept fixed. 

So far we have seen how the electron star emerges continuously from the Reissner- 
Nordstrom black hole solution as the temperature is lowered below a critical value. The 
electron star is characterised by the fact that a nonzero fraction of the charge is carried 
by a fluid outside of the black hole. These two observations suggest that the electron star 
birth is a second or higher order phase transition in the system. To make this claim precise, 
we must compute the behaviour of the free energy across the transition. The free energy 
density of the theory is most easily obtained from the thermodynamic relation 

n = M-p,Q-§T. (17) 
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Here s is the (rescaled) Bekenstein-Hawking entropy density given by 

8^ L 2 



I? 

K 



2 S S K 2 4 r 2 



(18) 



It was shown in [2] that (17) held for zero temperature electron stars by showing that the 



Lagrangian ([!]) was a total derivative on shell. In the presence of a nonzero temperature 
horizon, it is easy to check that there is an additional contribution to the on shell action at 
the horizon that contributes the necessary sT term. 

Figure [4] compares the free energy of two electron stars with the Reissner-Nordstrom 
black hole in the absence of a fermion fluid. The first important observation is that the 
free energy of the stars is indeed lower than that of the black hole for all temperature 
below the transition temperature at which the star is born. Secondly, we see that the 
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Figure 4: Free energy of the Reissner-Nordstrom black hole (top, red) and the free energy 
of two electron stars (blue) as a function of temperature. The electrons stars have (5 = 10 
and rh = 0.2 (lower) and rh = 0.3 (upper). The maximal temperature plotted is the critical 
temperature of the lower electron star. 



transition is extremely soft. In the plot one can only distinguish the free energy of the black 
hole with and without the fermion fluid at temperatures below around half the transition 
temperature. Increasing the accuracy of the numerics and fitting the curves close to the 
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transition point, we find that the transition is third order, so that AQ ~ (Tq — T) 3 . This 
fact can be understood analytically as follows. 

Consider a fermion with mass rh < 1. We noted in figure [I] that a fermionic fluid will 



start to condense at a radius r = r c and temperature T = Tq given by the conditions (11). 
Cooling below the critical temperature, we saw in figure [2] how the fermion fluid spreads 
out as the temperature is lowered. At any temperature the fluid has width Ar = r% — r\. 
Solving the conditions /iio C .( r i,2) = m perturbatively at small T — Tq, using the background 
Reissner-Nordstrom geometry, gives to leading order 

^ = #(l-T/T c ) 1 / 2 + ... . (19) 

Here and below $ refers to a complicated but computable mass dependent number. One 
can further check that to leading order in T — Tq the backreaction of the fermion fluid on 



the geometry neither alters nor contributes to the result (19). Because the backreaction of 
the fluid on the black hole solution appears to be subleading, we might suspect that the 
contribution of the fluid to the free energy is simply, from the Lagrangian ([!]), 

Ah = TAS E = - p ^Tglr 2 dr = -# /3 T 3 (1 - T/T c f . (20) 

J r\ 

Here Se is the rescaled Euclidean action, divided by the volume of the field theory spa- 



tial directions. The scaling follows immediately from (19) together with the facts that, 
from (jsj), we have p ~ $Pq/,c f° r Aloe. = rh + 5fi\ oc , and that <5/ii oc . = O ((Ar) 2 ). This 
second statement in turn follows from the fact that the local chemical potential may 
be approximated, again simply using the Reissner-Nordstrom geometry, by the parabola 
Aloe. = m + #(1 — T/Tc) 1 / 2 Ar/r + — #(Ar/r + ) 2 just below the transition temperature. To 
determine the coefficient in ( 20 ) , one should use the expression ^ for the pressure together 
with the parabolic expression for jii oc . and perform the integral. This can be done either 



numerically or analytically. Some illustrative values for the coefficient # in (20) are 

771 = 0.8 



rh = 0.1 rh = 0.3 
# = 0.0139 # = 0.7489 



rh = 0.5 
# = 3.529 



(21) 
# = 27.02 



We have checked that these values agree to within around one part in hundred with the 
values obtained by fitting the output of a full numerical integration of the equations of 
motion (using MATHEMATICAL NDSolve with WorkingPrecision set to 30). This provides 



good evidence that ( 20 ) is in fact the exact leading order result for the free energy difference 
close to the transition temperature and that there is therefore a third order phase transition 
in the system. 
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The reason the transition is so soft is a combination of two effects: immediately below 
the transition temperature there is a low density of fluid and it is furthermore constrained 
to a thin band of radii in space. While the local charge and energy densities both scale 

3/2 5/2 

like <5Aio C close to the critical temperature, the pressure scales like <5/t loc . This additional 
suppression feeds through to the fluid action and hence the free energy. 

Beyond the free energy, a second interesting thermodynamic observable is the entropy 
density. The original motivation to consider the backreaction of fermions on black hole 
spacetimes was to replace the emergent IR scaling of extremal black holes, which has an 
infinite dynamical critical exponent z = oo, with a finite z scaling geometry at zero tem- 
perature p]. Dimensional analysis implies that at low temperatures, T <C /t, the entropy 
density scales as s ~ T 2 I Z . This scaling has previously been exhibited in a different holo- 
graphic system in [16]. In [2] the exponent z was computed in terms of {/3,m}. Reading off 
an entropy scaling from our numerics therefore gives a vivid diagnostic of low temperature 
criticality as well as a check on our numerics. Figure [5] shows the entropy as a function 
of temperature for the black hole in the absence of an electron fluid together with three 
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Figure 5: Temperature dependence of the entropy density. Red line (top) is Reissner- 
Nordstrom. The three blue lines (lower) are electron stars with j3 = 20 and, from left to 
right, rh = 0.7, 0.36, 0.1. Fitting to a power law at low temperatures leads to z ~ 5.4, 2, 1.5 
respectively. 
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electron stars. The entropy density of the black hole tends to a constant at zero tempera- 
ture while the electron star entropy densities tend to zero. By fitting the numerical output 
to a power law, we can read off the dynamical critical exponent via s ~ T 2 I Z . The values 
obtained, quoted in the caption of figure [IJ agree with the zero temperature results of [2] 
to within one part in a hundred. 

Two issues remain to be discussed. Firstly, we promised a more detailed explanation of 
why it is consistent to treat a zero temperature fermion fluid in a black hole background. 
The black hole implies that the Euclidean time circle is periodically identified with radius 
1 /T. The local charge density in ^ should therefore necessarily include thermal excitations 

2 f d 3 p ±1 

a = 2tt fi^J (27r)3 - - e(B z FWoc .)/T loc . • {22 > 



The local chemical potential is again (|4|). The local temperature T\ oc = T '/ \fgu = T/Ly/J. 
It is then clear from the definitions and our limit ^ that 

^ = 0f4)«l. (23) 
Woe. \eLJ 

This justifies treating the local fermion fluid at zero temperature as we have done through- 
out. The physical content of this approximation is that while fermions are inherently quan- 
tum mechanical, the large local chemical potential allows the fermions to be treated in 
terms of local classical fields, and therefore suppresses the effect of Hawking radiation. As 



usual, the relation (23) will break down parametrically close to the black hole horizon due 
to the diverging local temperature. A proper treatment of the quantum stress tensor will 
indicate that there is no physical divergence at the horizon [17j . 

Finally, while the low temperature phase can be characterised in the bulk by the fact 
that a nonzero fraction of the total charge is carried by a fermion fluid, as seen in figure [3| 
we have not discussed an order parameter of the boundary quantum field theory that could 
distinguish the two phases. The bulk microscopic fermion field is dual to a 'single trace' 
fermionic operator ^ in the dual quantum field theory. It is useful to Fourier decompose 
this operator into creation operators ■ [H |9l [3] . On shell bulk fermionic states have 
uj = u>fc(M), where M is a label that depends on the radial location of the bulk fermion 
as well as its momentum in the radial direction. Details are given in [3] where it is also 
(implicitly) noted that the charge carried by fermions can be expressed as the expectation 
value of the 'double trace' operator 

Qfermion. = / dMA(M) J Sk (<£ {M) ^ . (24) 
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The function A{M) is to be read off from the distribution of fermions in the bulk. This 
quantity appears not to be easy to express purely in terms of boundary quantities. A more 
natural boundary quantity is the generalized fermion density 

n = j ( c l,k c ^k) ■ (25) 

This expectation value is computed as follows. Each bulk fermionic state corresponds to 
a solution of the Dirac equation in the background electron star spacetime. Although 
within our Thomas-Fermi description the local fermion density drops to zero at the outer 
boundary of the electron star, each occupied state will have a nonvanishing tail that reaches 
the boundary of the spacetimej^] Squaring the coefficient of each {to, k} mode will determine 
the expectation values in ( |25[ ) via the usual holographic dictionary. It should be possible to 
perform this computation explicitly and we hope to return to this point in the future. 

To leading order in the bulk semiclassical limit Q, the generalized fermion density 
(25) acts as a field theory order parameter; it is nonzero in the low temperature phase 
but vanishes above Tq- Away from the strict semiclassical limit it is no longer an order 
parameter. Thermal excitations imply that the generalized density is never exactly zero. 
The 'order parameter' does not break any symmetries of the theory. A third order transition 
is somewhat exotic, and deserves a better field theoretic understanding. 
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Noted Added. Yesterday the paper [18] appeared which has a strong overlap with our 
results. While mostly in agreement, that paper claims the transition is second rather than 
third order. We believe this is because they have neglected a — sT term in the free energy. 
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